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ABSTRACT 

The central kiloparsec region of the Andromeda galaxy is relatively gas poor, while the interstellar medium appears 
to be concentrated in a ring-like structure at about lOkpc radius. The central gas depletion has been attributed to a 
possible head-on collision 200 Myr ago, supported hy the existence of an offset inner ring of warm dust. We present 
new IRAM 30 m radio telescope observations of the molecular gas in the central region, and the detection of CO and its 
isotopes '^CO(2 - 1) and C**0(2-l), together with the dense gas tracers, HCN(l-O) and HCO-t-(l-O). A systematic 
study of the observed peak temperatures with non-LTE equilibrium simulations shows that the detected lines trace 
dense regions with in the range 2.5 x 10"* - 5.6 x 10^ cm“^, while the gas is very clumpy with a beam filling factor 
of 0.5 - 2 X 10“^. This is compatible with the dust mass derived from the far-infrared emission, assuming a dust-to-gas 
mass ratio of 0.01 with a typical clump size of 2 pc. We also show that the gas is optically thin in all lines except for 
'^CO(l-O) and *^CO(2-l), CO lines are close to their thermal equilihrium condition at 17-20 K, the molecular hydrogen 
density is larger than critical, and HCN and HCO-l- lines have a subthermal excitation temperature of 9 K with a density 
smaller than critical. The average '^CO/'^CO line ratio is high (~ 21), and close to the '^CO/C'*0 ratio (~ 30) that was 
measured in the north-western region and estimated in the south-east stacking. The fact that the optically thin *^CO and 
C**0 lines have comparable intensities means that the secondary element '^C is depleted with respect to the primary 
'^C, as is expected just after a recent star formation. This suggests that there has been a recent starburst in the central 
region, supporting the head-on collision scenario. 

Key words, galaxies: abundances; galaxies: individual: M31; galaxies: kinematics and dynamics; submillimeter: ISM; 
molecular data 


1. Introduction 

While the outer parts of the Andromeda galaxy exhibit 
numerous relics of past interactions (McConnachie et al. 

2009) , the main disc exhibits an enhanced star formation 
activity in a ring at lOkpc (Ford et al. 2013) and the cen¬ 
tral part is atypical. The M31 galaxy hosts a very mas¬ 
sive black hole with a mass of 0.7 - 1.4 x 10 *^Mq (Bacon 
et al. 2001; Bender et al. 2005) that is one of the most 
under-luminous supermassive black holes (Garcia et al. 

2010) , which is surrounded by very little gas (Melchior 
& Combes 2013). Infrared data and optical ionised gas 
display a 0.7 kpc off-centred inner ring, and Block et al. 
(2006) have proposed a frontal collision with M32 to ex¬ 
plain this double ring structure. This collision would have 
occurred 200 Myr ago. In this paper, we present observa- 

* Based on observations carried out with the IRAM 30m radio 
telescope. IRAM is supported by INSU/CNRS (France), MPG 
(Germany) and IGN (Spain). 


tions of millimetre molecular lines probing dense gas re¬ 
gions along the minor axis in the inner ring. In Sect. 2, we 
present the data used in this paper. In Sect. 3, we describe 
the data reduction. In Sect. 4, we describe our results and 
which physical information we can extract. In Sect. 5, we 
argue that our results support a 200 Myr old starburst trig¬ 
gered by the frontal collision with M32. For consistency 
with our previous analysis, we consider a distance of M31 
of 780kpc and a pixel size 1 arcsec = 3.8 pc. 

2. The data 

2.1. Millimetre observations 

We observed five positions listed in Table 1, which also 
gives the offsets and integration time. We computed the 
offsets with respect to the position of the optical cen¬ 
tre ((J2000) 00H2'”44.37F h- 41°16'08.34") provided by 
Crane et al. (1992) and used in previous papers (Melchior 
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Fig. 1. Positions observed in HCO+(1-0), HCN(l-O) and '^CO(2-l) (and C**0) are indicated with circles superimposed on the 
PACS/Flerschel lOOyum map (Viaene et al. 2014). Contour levels (0.05, 0.1, 0.2) correspond to Ag extinction computed assuming 
all the dust is in front of the bulge as discussed in Melchior et al. (2000). Extinction is underestimated by an unknown factor due to 
the uncertainties in the geometrical configuration of the gas in this area. 


Table 1. M31 observed positions. 


Pos 

RA DEC 

(J2000) (J2000) 

Aa 

(”) 

A(5 

(”) 

T/or 

(min) 

I 

00 42 29.1 -h41 18 03.6 

-172.1 

115.3 

284 

I-B 

00 42 26.4 +41 17 58.0 

-202.5 

109.7 

173 

2d-a 

00 42 56.2+41 14 38.0 

133.4 

-90.3 

262 

2d-b 

00 42 57.3 +41 14 38.0 

145.8 

-90.3 

204 

2c-a 

00 42 54.9 +41 14 14.0 

118.8 

-114.3 

79 


Notes. We provide the J2000 coordinates, offsets and integration 
times of the five positions observed. 


et al. 2000; Melchior & Combes 2011, 2013). The ob¬ 
served positions lie in the inner 1 kpc ring along the mi¬ 
nor axis (of the main disc inclined at 77 deg). The total 
integration time for each position is provided in the last 
column. In Figure 1, these positions are superimposed 
on the Herschel 100/tm map provided by Viaene et al. 


(2014) and extinction contours computed in Melchior 
et al. (2000). In the 100/tm and As maps, we can see the 
0.7 kpc off-centred inner ring. 

The observations have been performed with the 
IRAM 30 m radio telescope equipped with the EMIR re¬ 
ceiver in the period 20'^'-23''‘^ December 2011. They were 
performed in dual polarisation and the E0/E2 band com¬ 
bination were tuned in order to cover the following fre¬ 
quency band widths: 

1. 87.4-91.4GHz and 83.8-87.8GHz (optimised for 
HCO-h(I-O) and HCN(l-O)) 

2. 218.65-222.65 GHz and 214.95-218.95 GHz (opti¬ 
mised for 13 cO( 2-1) and C**0(2-l)). 

Two backends were connected to two main receivers: ETS 
with a spectral resolution of 0.19MHz and WILMA auto¬ 
correlator with a spectral resolution of 2MHz. WILMA 
was non-operational in some observations, and have 
mainly been used to check the quality of the observations. 
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Table 2. Dust properties derived from the infrared data for the M31 observed positions. 



D14 

V14 

Pos 

tJsPIRE350 

Um160 

'ycoW dust 

^SPIRE350 

(K) 

'T'CO/rf dust 

^M160 

(K) 

Tt 

dust 

(K) 

rpVV 

dust 

(K) 

(10^ Mo) 

M. 

(10’Mo) 

SFR 

(Mo/Myr) 

Tv 

I 

1.53±0.21 

2.1±0.5 

19.3 

20.3 

20. 

42. 

2.90 

3.2 

3.3 

0.52 

I-B 

1.52±0.02 

1.3±0.3 

19.3 

18.7 

20. 

57. 

2.86 

2.5 

2.8 

- 

2d-a 

2.73±0.37 

2.6±0.5 

21.3 

21.1 

23. 

34. 

0.98 

3.2 

13.4 

_ 

2d-b 

1.78±0.15 

2.6±0.5 

19.8 

21.1 

23. 

34. 

0.98 

3.2 

13.4 

- 

2c-a 

2.07±0.40 

3.1±0.4 

20.3 

21.8 

23. 

- 

1.49 

3.5 

16.1 

1.89 


Notes. In the second and third columns, we provide the mean starlight heating rates provided by Draine et al. (2014, D14) for two 
different resolutions (24."9 and 39") at the pointing positions. In the fourth and fifth columns, we present cold dust temperatures 
derived from these mean starlight measurements. In the sixth through tenth columns, we provide cold and warm dust temperatures, 
dust and stellar masses, star formation rates and dust optical depths computed by Viaene et al. (2014, V14) at these positions in 
36"pixels. 


Rest frequency (MHz) 


89400 89200 89000 88800 88600 



Fig. 2. Averaging of the observations performed in the 88 GHz 
bandwidth. The main astronomical lines present in the band¬ 
width are indicated (in black) in the rest-frame frequency of 
the Andromeda galaxy, whose systemic velocity is around - 
300kms“'. For each line, the redshifted and blueshifted parts 
of the expected velocity range are indicated. The noise level 
achieved is 0.2 mK for a frequency binning Av = 14 MHz. There 
are clear detections of HCO-)-(l-O) and HCN(l-O), with possible 
components at both sides of the systemic velocity. 
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Fig. 3. Averaging of the observations performed in the 216 GHz 
bandwidth. The '^CO(2-l) and C'*0(2-l) lines present in this 
bandwidth are indicated in the rest-frame frequency of the 
Andromeda galaxy. For each line, the redshifted and blueshifted 
parts of the expected velocity range are indicated. The noise 
level achieved is 3mK for a frequency binning Ay = 7 MHz 
(Av = 9.5kms“*). Apart from the multiple-component *^CO(2- 
1) signal, the C'*0(2-l) signal is weak in the receding part of 
the spectra. 


We smoothed the FTS data to the optimised resolution of 
2.1 km s ' up to 34km s *. Observations were performed 
in non-symmetric wobbler mode to avoid sky subtraction 
from Andromeda’s main disc. The azimuthal radius of the 
wobbler was set to 240 arcsec or set manually, relying on 
the extinction map of the area (Melchior et al. 2000). As 
discussed by Greve et al. (1996), large wobbler throws 
reduce the sensitivity to point sources. However, this is 
restored to a large extent with an active near-focus correc¬ 
tion. It is also possible that some emission is subtracted by 
the OFF signal. In principle, this would occur at a velocity 
different from the signal velocity, even though this cannot 
be excluded given the complex velocity pattern. One can 
also argue, on the basis of the detections discussed in the 
next section, that the emitted signal is weak, so the prob¬ 
ability of chopping a strong source is small. As shown 
in the figures discussed below, there is no obvious sign 
of emission in the OFF measurements, so this is in the 
worst case a second order effect. In addition, these effects 


are expected to affect *^CO(2-l) and C'*0(2-l) similarly, 
as well as HCOh-(I-O) and HCN(l-O). They could never¬ 
theless affect the normalisation with respect to *^CO(2-l) 
and *^CO(1-0) obtained independently. 

We retrieved '^CO(2-l) and '^CO(l-O) data from 
previous IRAM observations described in Melchior & 
Combes (2011, 2013) (and in prep.). *^CO(1-0) detection 
is only available for M3 II. 

2.2. Infrared data 

The five positions observed at IRAM 30 m were observed 
by Herschel/SPIRE instrument. The 1-kpc ring hosts cold 
temperature gas, while inside the ring the gas is hotter and 
not detected at large wavelengths. We extracted different 
parameters relevant for our study from the work of Viaene 
et al. (2014) and Draine et al. (2014). These parameters 
are summarised in Table 2. 
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Fig. 4. Carbon monoxide spectra observed on the north-west 
side of the inner ring are superimposed in velocity. The main 
beam temperature is provided for ’^CO, while the *^CO and 
C'*0 spectra are multiplied by a factor of 10 and 4 in the bottom 
and top panels, respectively. 



Fig. 5. Carbon monoxide spectra observed on the south-east side 
of the inner ring are superimposed in velocity. The main beam 
temperature is provided for '^CO, while the '^CO and C'*0 
spectra are multiplied by a factor of 2, 10, and 1 for the top, 
middle, and bottom panels. 


Draine et al. (2014) adjusted a physical dust model 
on Herschel and Spitzer observations to estimate the 
mean intensities of starlight heating the dust. These au¬ 
thors computed these mean intensities for the SPIRE 
350/Herschel and MIPS 160/Spitzer resolutions of 24."9 
and 39", respectively. Viaene et al. (2014) account for 
two more photometric channels. The mean starlight re¬ 
covered for both resolutions are provided in the second 
and third columns of Table 2. These differences have an 
amplitude comparable to the 30% discrepancy between 
the inferred heating rate and the value predicted by a sim¬ 
ple model of bulge starlight discussed by Draine et al. 
(2014). This provides an idea of the real uncertainty of 
these estimates, which smooth the substructures, while 
the corresponding resolutions (24."9 and 39") are compa¬ 
rable and larger than those of our millimeter observations 
(11", 21", and 28"; see Table 3). We then provide the 
cold dust temperatures computed according to the Draine 
& Li (2007) model as In the studied area, 

Draine et al. (2014) find that the starlight heating inten¬ 
sity is compatible with the stars in the bulge, as found by 
Groves et al. (2012). Relying on Herschel data, as well as 
GALEX, SDSS, WISE, and Spitzer, Viaene et al. (2014) 
model the spectral energy distribution to extract physical 
quantities. They estimate cold dust temperatures, which 
are similar to those derived from Draine et al. (2014), as 
shown in Table 2. The differences could be accounted for 
by the different datasets used as well as different mod¬ 
elling, Viaene et al. (2014) directly fit the spectral en¬ 
ergy distribution, while Draine et al. (2014) relies on a 
physical model. Viaene et al. (2014) have also provided 
a warm dust temperature (labelled ’’birth cloud” in their 
paper) larger than in the surrounding regions. This sug¬ 
gests that there are substructures linked to recent star for¬ 
mation, corresponding to (weak) star-forming activity (3- 
15Mo/Myr). We can expect the kinetic temperatures to 
be equal at most to the dust temperature, depending on 
the molecular hydrogen density. In the following, we con¬ 
sider a kinetic temperature of 20 K, which corresponds to 
the cold gas temperature. 


3. Data reduction 

We used the CLASS' package for the data reduction. 
We inspected individual spectra. Various interference 
lines are present far from the expected astronomical sig¬ 
nals. These interference lines affect the determination of 
the baselines, and they have been carefully removed on 
each spectra. The M3 II position has suffered an inter¬ 
ference next to the strongest peak of the '^CO(l-O) line. 
This interference only affects the data acquired on 20'^’ 
December 2011 (18/30). We thus set to zero the corre¬ 
sponding channels on data taken with the ETS backend 
with the best spectral resolution (0.26kms^*). 

The observations averaged in the two main band- 
widths are presented in the Eigs. 2 and 3. We used default 

' Continuum and Line Analysis Single-dish Software, 
http://www.iram.fr/IRAMFR/GILDAS 
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Fig. 7. Stacking centred at -144.3 km s“’ of each *^CO(2-l) velocity component with v > -300kms“*. The '^CO component is 
displayed in red (at -144kms“'). There is a component (displayed in green at -1285 kms^*) in C'^O which appears at the same 
velocity. The lines adjustment are displayed in blue. The velocity ranges corresponding to the two lines are displayed as rectangles. 
Standard weights (oc Ti„ dvjTsys^) for each observation have been used for the averaging. 
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Fig. 6. Signal detected in HCO+(1-0) and HCN(l-O) at the po¬ 
sition of M3 II (top) and M31I-B (bottom). The main-beam tem¬ 
peratures are displayed as a function of the rest-frame frequency 
of the Andromeda galaxy (corresponding to -300 kms^*). 


weighting with exposure time and system temperature. 
These spectra reveal multicomponent signals correspond¬ 
ing to the Andromeda velocity range (given its systemic 
velocity of -300kms“'). 

Characteristics of the main astronomical molecu¬ 
lar lines (Schoier et al. 2005) present in the observed 
bandwidths are provided in Table 3 as follows: half 
power beam width (HPBW), beam efficiency, rest-frame 


frequency, and velocity bin corresponding to 14 MHz. 
We also provide the following characteristics of these 
molecules: energy level, Einstein coefficient Auk, colli¬ 
sional de-excitation rate and critical density n^ com¬ 
puted at 20 K. If we consider a collisional temperature of 
50 K, the critical densities of HCN and HCO-h increase 
by 56% and 21%, respectively, but by 7.4% for CO. As 
indicated in Figs. 2 and 3, there are obvious detections 
of HCOh-(I-O), HCN(l-O), and '3cO(2-l) in the summed 
spectra and a weak signal in C'*0(2-l). There is no ob¬ 
vious contamination from the Milky Way, which could be 
possible as M31 is 20 deg below the Galactic plane. The 
spectra obtained in the lower outer bands are presented in 
Fig. A.l but do not exhibit any signal. They are used to 
derived upper limits. For each position, the detections are 
adjusted with a Gaussian function. Up to three different 
CO lines are detected at different central velocity at each 
position. As discussed in Melchior & Combes (2011), this 
is due to several intervening molecular components along 
the line of sight. 

All of the measurements are summarised in Table 4. 
The columns on the right side provide the CO line mea¬ 
surements. We provide, for each CO line, the central FSR 
velocity, the FWHM velocity dispersion Av, the peak 
(main-beam) temperature 'Ipeak, and the integral I^o of 
the line. In the last two columns, the '^CO/'^CO(2-l) and 
'^C0/C^®0(2-1) line ratios and lower limits are shown. 
We also provide (in the bottom of this table) our '^CO(l- 
0) measurements described in Melchior & Combes (2011) 
for M31-I with a beam resolution of 21 arcsec, similar 
to the resolution achieved for HCN(l-O) and HCO-i-(l- 
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0), and the column density derived as discussed be¬ 
low but directly computed from '^CO(l-O). These column 
density estimates differ by up to a factor of 3 with respect 
to '^CO(2-l)-based estimates, due to the difference in res¬ 
olution and corresponding dilution effects. 

The columns on the left side of Table 4 list HCN(1- 
0) and HCOh-(I-O) measurements and 3cr upper limits 
(for Av=32kms“') for some of the lines presented in 
Table 3. The central table provides the remaining upper 
limits, as well as column densities (and upper limits) of 
the species with some detections, namely '^CO, C^^O, 
HCN, and HCO-h, assuming local thermal equilibrium 
(LTE) and optically thin conditions. The column densi¬ 
ties have been corrected for the beam filling factor, as 
further discussed in the next section. Given the veloc¬ 
ity dispersion of the observed HCN lines (>25kms '), 
the HCN hyperfine line structure (~ 10 km s ') will not 
be detected. We also provide the N^j column densities 
of molecular hydrogen relying on Xco and CO(2-l)/(l- 
0) line ratio values. We assume a Galactic ratio Xco - 
NH2IIC0 = 2. X 10^'’cm“^(Kkms“')“' following Strong 
et al. (1988) and Dame et al. (2001), which is compatible 
with the findings of Smith et al. (2012) for Andromeda 
and Mauersberger et al. (2003) for local galaxies. As fur¬ 
ther discussed in Sect. 4.3, we take = Ico/Vbf- We 
assume a CO(2-1)/(1-0) line ratio of 0.8, as measured by 
Melchior & Combes (2011) in a clump located within 
10"of M31I (and named M31G), to compute Ahj column 
densities. 

4. Results 

In this section, we describe our results. In Sect. 4.1, we 
describe the detections of five different molecular emis¬ 
sion lines in the two observed positions along the mi¬ 
nor axis. In Sect 4.2, we compute the line ratios relying 
on previous detections and compare them with previous 
results achieved in other galaxies, in particular. Galactic 
giant clouds and in M3rs disc. In Sect. 4.3, we discuss 
the basic properties we can derive from the detections, in¬ 
cluding excitation temperatures, optical depths, beam fill¬ 
ing factors, and column densities, and the uncertainties 
and limitations of this type of analysis. In Sect. 4.4, we 
discuss the parameters estimated under LTE conditions 
and compare them with Galactic values. In Sect. 4.5, we 
run RADEX (van der Tak et al. 2007) simulations to relax 
the LTE conditions and test the validity of our previous 
assumptions. While CO emission lines correspond to gas 
in thermal equilibrium, gas traced by HCN and HCO-i- 
emission lines is under subthermal excitation conditions. 

4.1. Description of the detections 

The CO lines detected in the north-west (south-east) re¬ 
gion are shown in Eig. 4 (Eig. 5). The HCN(l-O) and 
HCO-i-(l-O) lines detected on the north-west side are 
shown in Eig. 6. As shown in Table 4, we detected up 
to five different molecular lines (two to four per posi¬ 
tion and per velocity component). Eor some positions 


(M3 II, M31I-B, M31-2d-a), there are clearly several 
well-identified peaks with different velocities, while for 
other positions the signal is very wide and best fitted with 
several components. The most striking cases are M31-2c- 
a, and the low-velocity component of M31 -2d-b. The 4cr 
*^CO(2-l) signal measured for M31-2c-a is affected by a 
small integration time and the detection integrates several 
*^CO(2-1) components. Hence, this signal is not used for 
subsequent stacking. Apart from the *^CO and some *^CO 
lines, the detections are weak with peak temperatures at 
the mK level. We detected HCO-h(I-O) (HCN(l-O)) at the 
6cr and 9cr (4cr and 5cr) levels, while most CO(2-l) de¬ 
tections are above 4cr with a strong one at Her; but three 
weak *^CO detections in the south-east region are kept at 
the 3.3cr level. Knowing the low intensities of these lines, 
we can infer that it is extremely challenging to perform 
multi-transition observations of these molecules. 

While '^CO and '^CO are detected on both regions, 
C'^O, HCN(l-O) and HCO-i-(l-O) are only detected on the 
north-west side. If the dust-to-gas mass ratio is similar on 
both sides, according to Viaene et al. (2014) we expect a 
gas mass a factor of 3 larger on this side (see Table 2). 
Similarly, when one averages the observed column densi¬ 
ties of molecular hydrogen provided in Table 4 along the 
two lines of sight, the north-west positions have a column 
density a factor of 2 higher than the south-east positions. 
Also, the radiation field (U) is larger on the south-east side 
than in the north-west side: this may be related to larger 
CO column densities and to the fact that dense gas (C'^O, 
HCN, HCO-I-) is not directly detected in the south-east 
positions. 

The C'^O detection corresponds to the component 
at the systemic velocity, while the HCN and HCO-t- de¬ 
tections are in the blueshifted components. However, as 
shown in Eig. 7, stacking of the v > -300km s ' '^CO 
detections (excluding the M31-1(3) detection) also reveals 
a redshifted component in C'^O with an amplitude com¬ 
parable with '^CO. The characteristics of this stacking 
detection is presented in Table 5. We thus have at least 
two different C'^O components in two different velocity 
ranges. 

Eigure A.2 denotes the stacking of v< -300km/s 
and v> -300km/s of HCN(l-O) and HCO-h(I-O) spec¬ 
tra. The detected signals are summarised in Table 5. The 
south-east spectra have been shifted to the '^CO velocities 
and averaged with standard weights (oc Ti„, dvIT^ys^). Eor 
HCN, there is 3.7cr detection in the blue stacking while 
there is a possible 2.6cr detection in the red stacking. In 
contrast to the direct observations, the stacked HCN sig¬ 
nal is stronger than the HCO-h signal, with a HCO- 1 -/HCN 
line ratio of 0.87 + 0.49 in the red stacking and < 0.68 in 
the blue stacking. 

The difficulty of the interpretation of these dense gas 
detections is due to a complicated configuration. We inte¬ 
grate the signal over a relatively large beam (42-106 pc in 
projection), and the velocity dispersions are significantly 
larger than those (~2km/s) detected in M31’s main disc 
(Schruba et al., in prep) with Av > 30kms“' for most 
of our detections. This is because of the strong velocity 
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Table 3. Properties of the main lines in the two observed bands. 


Transition 

HPBW 

(arcsec) 

^mb 

Vo 

(GHz) 

^v(lMHz) 

(kms^') 

E„p 

(K) 


Cul 

(cm" s-') 

n, (20 K) 
(cm-") 

SiO 7 = 2 ^ 1 

28. 

81 

86.847 

3.5 

6.25 

2.93 10^" 

Tr™ 

2.93 10" 

C2HA= 1 -»0 
(7 = 3/2^ 1/2,F= 1-1) 

28. 

81 

87.329 

3.4 

4.19 

0.26 lO^*' 

5.1910-'" 

5.01 10“* 

HNCO(4o4 ^ 3o3) 

28. 

81 

87.925 

3.4 

10.55 

9.02 10-*’ 

9.210-'" 

9.80 10" 

HCN y = 1 0 

28. 

81 

88.632 

3.4 

4.25 

2.4110-^ 

1.9210-“ 

1.2610° 

HCO+ y = 1 0 

28. 

81 

89.189 

3.4 

4.28 

4.25 10-^ 

2.3 10-'° 

1.8510" 

HNC y = 1 0 

27. 

81 

90.664 

3.3 

4.35 

2.69 10-^ 

8.95 10-“ 

3.01 10" 

HCjNy = 10^9 

27. 

81 

90.979 

3.3 

24.01 

5.81 10-^ 

1.110-'° 

5.28 10" 

Sioy = 5 ^4 

11. 

60 

217.105 

1.4 

31.26 

5.20 10-^ 

1.110-'° 

4.73 10° 

H2CO(3o3 ^ 2o2) 

11. 

60 

218.222 

1.4 

21.0 

2.82 10-“* 

1.110-'° 

2.56 10° 

H2C0(322 ^ 22,) 

11. 

60 

218.476 

1.4 

68.1 

1.5710-“' 

5.3 10-“ 

2.96 10° 

H2C0(321 ^ 22o) 

11. 

60 

218.760 

1.4 

68.1 

1.58 10-“* 

4.3 10-“ 

3.67 10° 

c'^'cy = 2 ^ 1 

11. 

60 

219.560 

1.4 

15.81 

6.0110-^ 

6.44 10-“ 

9.33 10" 

‘^coy = 2 ^ 1 

11. 

60 

220.399 

1.4 

15.87 

6.04 10-^ 

6.44 10-“ 

9.38 10" 

‘^coy = 2^1 

11. 

60 

230.538 

1.3 

16.60 

6.9110-'' 

6.44 10-“ 

1.07 10“* 

‘^coy = 1 0 

21. 

71 

115.271 

2.6 

5.53 

7.20 10-" 

3.25 10-“ 

2.22 10" 


Notes. The beam sizes (HPBW) and the main beam efficiencies are presented with the rest frequencies Vo of the brightest 
lines in each band and the velocity bin corresponding to 1 MHz. The energy levels E„p and Einstein coefficients A„; are based on 
the Leiden Atomic and Molecular Database (Schdier et al. 2005); the collisional de-excitation rate coefficients C„( are computed at 
a kinetic temperature of 20 K. The critical density is computed as n^. = A„;/C„;. 


gradient present in this region and to the superposition of 
several clouds. This also points towards inner elongated 
orbits observed with a large radial velocity range (in con¬ 
trast to the outer orbits observed in the main disc where 
small velocity dispersions (~2kms“^) are expected and 
observed). M31-I(l) has a smaller velocity dispersion, 
which is still four times larger than disc clouds. The other 
components of M31 -I are very large and could correspond 
to an underlying (possibly continuous) component with a 
150 km s ' width with non-circular orbits. This could be 
a relic of a spiral arm or dissolved bar. 

4.2. Observed line ratios 

In this section, we analyse the observed line ratios, which 
are directly derived from the observations. 

When one computes line ratios based on velocity- 
integrated intensities as shown in Table 6, dilution effects 
due to the unknown beam filling factors should be re¬ 
moved at least at first order. Empirical comparisons with 
similar targets can be performed. 

Figure 8 shows that our measured *^C0/C**0(2-1) 
line ratios lie below the (1-0) measurements of Liu et al. 
(2014) in infrared Galactic dark clouds (average 4.4+1.3; 
see their Fig. 2), but corresponds to the lower end of 
the Tan et al. (2011) starburst sample distribution based 
on CO/C *^0(1-0) line ratio. Our measured line ratios 
are also significantly lower than the M51’s spiral arms 
gas detected by Schinnerer et al. (2010) (their average 
'3CO/C'*^0 line ratio: 3.5). 

It is remarkable that the C'^O detection is at a level 
comparable to the *^CO detection, in contrast to what 
is observed in the Galaxy (e.g. Nishimura et al. 2015; 


Shimajiri et al. 2014), where the C'^O is usually 10 to 
150 times weaker than CO. *^CO/C'*0 line ratios close 
to one are expected if *^CO is optically thick or if '^C is 
depleted with respect to '^C, such as after a starburst. 

In Fig. 9, we compare our *^CO/'^CO(2-l) line ra¬ 
tios with values measured previously in other galaxies. 
We consider the following catalogues: 

- the (2-1) line ratio compilations for different morpho¬ 
logical types from Krips et al. (2010) (We provide 
the different morphological types under the form of 
hatched histogram); 

- the (1-0) line ratio measured by Li et al. (2015) in 
nearby galaxies; 

- the (1-0) line ratio measured by Matsushita et al. 
(2010) in starburst and post-starburst galaxies; and 

- the (1-0) line ratio measured by Tan et al. (2011) in 
nearby infrared-bright galaxies and starbursts. 

The (1-0) line ratios have been corrected with the aver¬ 
age '2 cO/' 3CO(1-0)/(2-1) ratio (1.33) from Krips et al. 
(2010). These authors found that the '^CO/'^CO(2-l) ra¬ 
tios of three SAURON galaxies are somewhat higher than 
those in galaxies of different Hubble types. The hatched 
histogram of Fig. 9 shows the different morphological 
types when available, and indicates that there is no trend 
with morphology. As discussed by Li et al. (2015), the 
isotopologue ratio does not depend on the morphological 
type of the galaxy. Our measurements (with an average 
'^CO/'^CO= 20), shown with a thick line hatched his¬ 
togram, lie below the standard isotopic ratio (53) (Wilson 
& Rood 1994) but fall in the same range as the other mea¬ 
surements. This is consistent with the fact that *^CO is 
optically thick. As a result of very different cloud ge¬ 
ometries, it is difficult to compare our *^CO/*^CO(2-l) 
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Table 4. Results. 
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Table 5. Line measurements for the stacking detection shown in Fig. 7, based on the red (v> -300km/s) velocities of both regions 
(excluding the M31-I (3) detection and M31-2c-a(l) tentative *^CO(2-l) detection). 


Pos. 

'^S tack 

(Vo 

) (km s *) 

Av (km s ‘) 

Tpm* 

(mK) 

I/me (K km s') 



'^CO 

c"*o 

"CO 

C"0 

"CO 

C"0 

'^CO 

C'O 



(2-1) 

(2-1) 

(2-1) 

(2-1) 

(2-1) 

(2-1) 

(2-1) 

(2-1) 

Stack (red) 

-144.2 

-212 

-212 

50* 

50* 

1.65 

0.95 

0.088±0.015 

0.051±0.015 





31 




0.15±0.03 


SE stack (red) 

-226. 

-218. 

-218. 

39* 

39* 

2.26 

1.28 

0.094±0.015 

0.053±0.015 





39 




0.14±0.04 



Pos. 

'^S tack 

(Vo 

) (km s ‘) 

Av (km s ‘) 

Tpm* 

(mK) 

I;,ne(Kkms ') 



HCO+ HCN 

HCO-t- 

HCN 

HCO-t 

HCN 

HCO-h 

HCN 



(1-0) 

(1-0) 

(1-0) 

(1-0) 

(1-0) 

(1-0) 

(1-0) 

(1-0) 

SE Stack (blue) 

-523 

- 

-524.5 

89* 

89±25 

- 

1.17 

< 0.075 

0.11 ±0.03 

SE Stack (red) 

-234 

-231 

-231 

78±29 

78* 

0.74 

0.85 

0.061±0.026 

0.070±0.024 


Notes. We present the central velocities, velocity dispersions, peak temperatures, and integrated lines. (Vo) is the averaged velocity 
of the spectra used for each stacking. Similar stacking based on the '^CO(2-l) velocities is performed for FICN and FtCO+ spectra 
of the south-east (SE) positions. As shown in Fig. A.2, redshifted (blueshifted) velocities are considered separately. 


Table 6. Observed line ratios. 


Position 

Vo (km/s) 

HCO±/HCN 

(1-0) 

HNC/HCN 

(1-0) 

HOC+/HCO± 

(1-0) 

'^CO/HCO± 

(1-0) 

‘^CO/HCN 

(1-0) 

'^corco 

(2-1) 

1(2) 

~ -402 

1.5±0.5 

<1.09 

<0.59 

43.±8 

62.±18 

24.±5 

I-B(l) 

~ -470 

1.7±0.4 

<0.35 

<0.21 

J3.± 7 

89.± 18 

30.±3 


Position 

Vo (km/s) 

HCO±/HCN 

‘^CO/C'O 

‘"co/c"'o 



"CO/'"CO 



(1-0) 

(2-1) 

(2-1) 



(2-1) 

1(3) 

~ -291 


28.6±9.5 

1.48±0.51 



19.4±5.0 

Stack (red) 



30±12 

1.73±0.59 



17.4±4.1 

SE stack (red) 


0.87±0.49 

33±14 

1.77±0.58 



18.7±5.1 


Notes. These ratio are computed with the velocity-integrated intensities measured within similar beams for the same component. 
Numbers in italics correspond to tentative estimates: *^CO(1-0) are derived from '^CO(2-l) with a 0.8 line ratio, stacked '^CO(2-l) 
and '^CO(2-l) are computed as the sum of the detected lines; and '^CO(2-l) and C'*0(2-l) are measured on stacked spectra. 


line ratio measurements in M31 with studies performed 
in Galactic regions. For instance, Orion regions observed 
by Nishimura et al. (2015) or Planck cold dust clumps 
(Wu et al. 2012) exhibits relatively low '^CO/'^CO line 
ratio (~ 2 - 3) because the dense clouds fill the beam and 
'^CO is very optically thick. 

Figure 10 shows the HCN/HCO-h( 1-0) and 
'^CO/'^CO(2-l) ratios we measured in the compo¬ 
nents approaching us in the north-western positions and 
the south-east stacking detection, which is redshifted. 
We compare these line ratios with previous observations 
from the literature. By superposing our measurements 
on the wide range of galaxy types observed by Krips 
et al. (2010), we find that our detections seem to lie 
in the peculiar/merger region, rather than in the star¬ 
forming (M33) and starburst part of the diagram. The 
HCN/HCOh- line ratio is compatible with the average 
ratio computed for M31’s disc by Brouillet et al. (2005). 
Again, the Galactic clouds of Nishimura et al. (2015) lie 
on the right hand side of this figure and are difficult to 
compare with M31 measurements. Our measurements 
are not compatible with Seyfert galaxies, which makes 
sense because there is no sign of active galactic nuclei 


(AGN) in M31 (Li et al. 2011). Interestingly, on the 
one hand, the HCN/HCO-h line ratio detected here is 
compatible with the value (1.5) measured by Liu et al. 
(2013) in infrared dark clouds. The red stacking ratio 
has a larger value but it is affected by a large error 
bar. On the other hand, HCN/HNC(l-0) line ratio is 
systematically larger than that measured in dark clouds 
by these authors. The HCN(l-O) line detections enable us 
to estimate lower limits on the HCN/HNC(l-0) line ratio 
towards the positions in the north-western area (M3 II 
and M31I-B), which are 0.9 and 2.9, respectively. Figure 
A.3 compares the HCOh-/CO( 1-0) and HCN/CO(1-0) 
line ratios detected on the north-western side of the inner 
ring with other positions in M31’s main disc (Brouillet 
et al. 2005) and with other galaxies. Our detections are 
compatible with the Brouillet et al. (2005) measurements 
performed in M31’s disc, but lie in the lower end of the 
overall distribution of Krips et al. (2010). Figure A.4 
shows the detections and the upper limits of our observed 
CO line ratios. 

Our detections differ from measurements performed 
in M33’s giant molecular clouds (Buchbender et al. 
2013). They are compatible with line ratios available for 
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Fig. 8. Comparison of our '^C0/C'*0(2-l) line ratio measurements (M31-1(3) and stacking) corresponding to different clouds (i.e. 
with different velocities) with previous measurements in nearby galaxies from Tan et al. (2011). The hatched histogram displayed 
the original '^CO/C'*O(l-0) line ratio. 


clouds in MSl’s disc studied by Brouillet et al. (2005), 
but there is no information on the CO isotopic ratio in the 
disc. The main trends regard the '^C0/C**0 line ratio, 
which is depleted, while the HCN/HCO-i-(l-0) line ratio 
is very typical of what is observed in other local galax¬ 
ies .Liu et al. (2013) have studied infrared dark clouds and 
argue that the abundance ratio of HCN and HCO-h may be 
unaffected by the environment. Even though the 100/tm 
and maps shown in Figure 1 do not exhibit strong fea¬ 
tures, we argue in Sect. 4.3 that our detections correspond 
to dense cores. These ratios depend non-linearly on phys¬ 
ical parameters, such as optical depth, excitation temper¬ 
ature, density, and some lines are probably optically thick 
and not necessarily in thermal equilibrium. This is ex¬ 
plored further in Sects. 4.3, 4.4, and 4.5. 


4.3. Estimate of the beam filling factors 

The temperatures we have detected for '^CO(2-l) are 
well below the cosmological background, even though the 
molecular lines are detected in emission. This is due to 
the low surface filling factor of the gas detected in the 
beam. In this section, we try to derive estimates of the 
beam filling factors affecting our observations. These es¬ 
timates affect all absolute parameters (temperature, den¬ 
sity, abundance, column density, etc.) that could be de¬ 
termined later, and thus the possible interpretation of the 
state of the detected gas. 

As the average '^CO/'^CO line ratio (20) is well be¬ 
low the standard Galactic abundance ratio (53) of 

Wilson & Rood (1994) and the C'^O line intensity com¬ 
parable to the '^CO intensity, we consider that the '^CO 


line is optically thick. As the gas is detected in emission, 
we can then derive a minimal excitation temperature of 
about 3K from the intrinsic brightness temperature Tg, 
with the usual formula 


TB^(f{Te.)-fiTbg))x(l-e-^) ( 1 ) 


where 


hv 


f(T)^ 


exp{!^) - 1 


( 2 ) 


Relying on the modelling of infrared measurements, 
Viaene et al. (2014) and Draine et al. (2014) estimated 
a cold dust temperature of the order of 20 K (see Table 2). 
Since we measured the line ratio CO(2-1)/CO(1-0)~0.8 
(Melchior & Combes 2011), which is relatively high, we 
assume that the line is thermalised and that we can define 
a single excitation temperature. This makes sense for cold 
gas in the central part of a large galaxy. Most importantly, 
one has to take into account the beam filling factor rjhf, 
which is linked to the intrinsic brightness temperature as 
follows: 


Tb - Tpeakhljf (3) 

where Tpeak is the peak temperature measured for our de¬ 
tections, i.e. the observed diluted main beam temperature. 
Given our assumptions about the excitation temperature 
and the '^CO optical depth, we can derive rjbf from our 
'^CO(2-l) detections as 

r]bf — EpgQk! [f{Tg,)-f(Tbs)\ (4) 

If the molecular lines are close to thermal equilibrium 
around 20 K, the beam filling factor is below one percent. 
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Fig. 9. Comparison of our '^CO/''^CO line ratio with previous measurements in nearby galaxies. The hatched histogram displayed 
the '^CO/*^CO line ratios with available morphological types. The Li et al. (2011), Matsushita et al. (2010) and Tan et al. (2011) 
samples measured the (1-0) transition and have been corrected by a factor 1.33 as discussed in the text. The Krips et al. (2010) 
sample and our measurements were directly measured the (2-1) transition. The black full line histogram gathers all the four samples 
corrected to (2-1). The black thick line histogram in the bottom corresponds to the measurements performed in this paper. The 
average *^CO/'^CO detected in this paper tends to be a bit larger than the main distribution from other areas and other galaxies (but 
well below the standard isotopic ratio of 53 (Wilson & Rood 1994)). The largest value above 50 might corresponds to an optically 
thin cloud, but it is affected by a large error bar (30%, see Table 4). 


This order of magnitude is consistent with several differ¬ 
ent qualitative arguments, as follows: (1) The beam size 
is 44pc for '^CO(2-l), while typical interstellar filaments 
have a size of 0.1 pc (Andre et al. 2014). (2) We are de¬ 
tecting gas complexes, which are probably extended and 
patchy according to the optical image. (3) The large ve¬ 


locity dispersions also suggest that we are integrating gas 
along a large depth. (4) Last, if one considers a standard 
gas-to-dust relation N// = 0.94 x 10^' x Ay cm“^ (Bohlin 
et al. 1978), we estimate Ay = 2.6 and 1.4 for M31-1(2) 
and M31-IB(1) (with no beam-filling factor corrections), 
while lower limits on extinction estimated on the optical 
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Fig. 10. Comparison with other catalogues (1). HCN/HCO+(1-0) line ratio versus '^CO/'^CO(2-l) line ratio. The two north-west 
detections presented in this paper (M3II and M31I-B) together with the red south-east stacking measurement are superimposed 
on previous measurements, namely: a compilation of various types of galaxies from Krips et al. (2010), early-type galaxies from 
Crocker et al. (2012) and M33’s giant molecular clouds from Buchbender et al. (2013). The hatched area corresponds to the mean 
HCN/HCO-l-(l-0) line ratio obtained by Brouillet et al. (2005) in M3rs disc. Surprisingly, our detections, compatible with the M31 
disc measurements of Brouillet et al. (2005), lie in ’’the peculiar/merger region”. 


images are of the order of Ab - 0.2 (Melchior & Combes 
2011). This order of magnitude of the beam filling fac¬ 
tor is also compatible with values previously discussed 
in the literature, for example Carlstrom (1988) estimated 
a volume filling factor smaller than 10~^ in M82, while 
Berkhuijsen & Fletcher (2015) estimate a filling factor of 
approximately 4% in M31 and M51. 

In the following, we consider that ^^CO(2-l) is opti¬ 
cally thick to determine the beam filling factors at ther¬ 


mal equilibrium, as provided in Table 7. The values 
we estimate for rji,/ are consistent with Ginsburg et al. 
(2011), who discuss that giant molecular clouds consist of 
very low volume-filling factor (~ 5 x 10“^) high-density 
(n(H 2 )~ 3 X 10^ cm“^) clumps, as well as typical starburst 
galaxies (with //*/ < 0.01). 
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Table 7. Estimate the beam filling factor r]i,f(Tg^) for different excitation temperatures, assuming that the '^CO(2-l) line is optically 
thick. 


Positions 

rti,f(25K) 

ribf{2bK) 

W(^5K) 

776/(10 ff) 

W(5K) 

"TlSCO 

■rci*o 

'rHCO+ 

Thcn 

^MR 

^OC 

1(1) 

2.1x10^^ 

2.9 X 10-3 

4.2x10-3 

8.0x10-3 

3.6x10-' 

0.049 




2.57 

0.97 

1(2) 

7.4x10-3 

9.9 X 10-3 

1.5x10-3 

2.8x10-3 

1.2x10-' 

0.033 


0.021 

0.012 

1.77 

0.67 

1(3) 

3.9x10-3 

5.2 X 10-3 

7.7x10-3 

1.5x10-3 

6.6x10-3 

0.056 

0.038 



2.97 

1.12 

I-B(l) 

1.8x10-3 

2.4 X 10-3 

3.6x10-3 

6.7x10-3 

3.0x10-' 

0.033 


0.016 

0.016 

1.74 

0.66 

I-B(2) 

2.5x10-3 

3.3 X 10-3 

4.9x10-3 

9.3x10-3 

4.2x10-3 

0.089 




4.73 

1.78 

2d-a(l) 

5.1x10-3 

6.7 X 10-3 

1.0x10-3 

1.9x10-3 

8.5x10-3 

0.059 




3.12 

1.18 

2d-a(2) 

5.3x10-3 

7.1 X 10-3 

1.1x10-3 

2.0x10-3 

9.0x10-3 

0.051 




2.69 

1.02 

2d-b(l) 

8.9x10-3 

1.2 X 10-3 

1.8x10-3 

3.3x10-3 

1.5x10-' 

0.018 




0.97 

0.37 

2d-b(2) 

2.2x10-3 

3.0 X 10-3 

4.4x10-3 

8.3x10-3 

3.8x10-3 

0.063 




3.36 

1.27 

2d-b(3) 

1.9x10-3 

2.6 X 10-3 

3.8x10-3 

7.2x10-3 

3.2x10-3 

0.135 




7.16 

2.70 

2c-a(l) 

3.1x10-3 

4.1 X 10-3 

6.1x10-3 

1.2x10-3 

5.2x10-3 

0.071 




3.78 

1.43 

2c-a(2) 

2.7x10-3 

3.7 X 10-3 

5.4x10-3 

1.0x10-3 

4.6x10-3 







2c-a(3) 

5.0x10-3 

6.6 X 10-3 

9.8x10-3 

1.9x10-3 

8.4x10-3 








Notes. The next four columns provide the optical depths for '^CO(2-l), C'*0(2-l), HCO+(1-0), and HCN(l-O) assuming an exci¬ 
tation temperature of 20 K and a beam filling factor equal to that derived from *^CO, i.e. 776/(20 K). The next two columns provide 
the '^CO(2-l) column densities computed for two standard '^CO/'^CO line ratios for the molecular ring (MR) and Galactic centre 
(GC) from Wilson & Rood (1994). 


For a given excitation temperature, we can derive the 
optical depth corresponding to each optically thin line: 


give an initial idea of the physical configurations of the 
observed complexes. 


"^mol — 



Tpeak 1 


( 5 ) 


We can then compute the column density associated to 
each molecular line assuming an excitation temperature 
T ■ 

, 6 ) 

hc^ g 1 — e r/hf 

where g is the level degeneracy g - 5 for CO(2-l) 
molecules and g = 3 for HCN(l-O) and HCOh-(I-O), Q(T) 
is the rotational partition function and is the energy 
level provided in Table 3. We also consider the beam fill¬ 
ing factors defined in Table 7. The factor enables 

to account for non-optically thin conditions (e.g. Mangum 
& Shirley 2015). However, in practice, we have to assume 
an unknown abundance ratio to determine Tmoi if it is 
optically thick, which limits the use of corrected column 
density in this analysis. 

The optical depths thus computed and provided in 
Table 7 show that our hypotheses are consistent; '^CO 
is optically thin. In Sect. 4.4, we discuss abundances esti¬ 
mated assuming LTE conditions. In Sect. 4.5, we explore 
how to derive molecular hydrogen densities for lower ex¬ 
citation temperatures with RADEX simulations (van der 
Tak et al. 2007). 


4.4. Local thermal equilibrium conditions 

We consider here that the different detected lines corre¬ 
spond to some gas in LTE with an excitation temperature 
of 20 K. This enables us to derive physical parameters 
from the observations with an analytical formula and to 


4.4.1. Optical depths and column densities 

With these assumptions, we derived optical depths for 
'^CO(2-l), Ci*^0(2-1), HCOh-(I-O) and HCN(l-O) with 
the filling factor 776/(20 X), computed in Sect. 4.3. These 
lines are optically thin, as shown in Table 7. 

If we consider standard '^CO/'^CO line ratios from 
Wilson & Rood (1994) of 53 for the molecular ring (20 for 
the Galactic centre) and the optical depths provided 
in Table 7, we estimate an average optical depth for '^CO 
of 3.2 (1.2). The '^CO(2- 1) transition is indeed optically 
thick. These different optical depths discussed above are 
provided for each line of sight in Table 7. 

In the middle panel of Table 4, we provide column 
densities of our main detections (*^CO, C^^O, HCN, 
HCO-h), some upper limits for optically thin gas (t << 1) 
in LTE, and the column density of molecular hydrogen 
derived empirically from the '^CO line intensity with a 
Xco factor as discussed in Sect. 3. All column densities 
account for the beam filling factor rji,/. 

4.4.2. Abundance ratios 

In Table 8, we tentatively compute the abundances (i.e. 
the ratio of the column density of the molecule with re¬ 
spect to the column density of molecular hydrogen) for 
the direct and stacking detections. In order to calculate 
the ratio of lines measured with similar HPBW resolu¬ 
tions, we use '^CO(2-l) to derive the molecular hydrogen 
column density for the '^CO(2-l) and C'^0(2-l) measure¬ 
ments and '2CO(1-0) (available for M31-I) for HCN(l-O) 
and HCO-i-(l-O) abundances. We thus study M31I(3) and 
M311(2). We also include the red stacking detections for 
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Table 8. Abundance ratios. 



Measured values (LTE) 

Galactic values 


1(3) 

1 Red Stacking 


Transition 

(2-1) 

/Mol 

X^^ '^CO 

/Mol 

(MR) (GC) 

‘“CO 

/Mol 

(MR) (GC) 


10.8+2.7x10-“ i 

7.3+2.4x10-“ 1.48+0.51 

&0±i.5x 10-“ i 

4.(5±i.4xl0-“ 1.73+0.58 

5.3x10-' 1.4x10-° 

8.6x10-“ 1.1x10-’ 

1 1 

6.2 12.7 


Measured values (LTE) 

Galactic values 


1(2) 

1 SE red Stacking 


Transition 

(1-0) 

Moi 

/HCO+ 

X^^ Moi 

/HCO+ 

Moi 

/HCO+ 

HCN 

HCO+ 

1.8+0.3 10-“’ 1.11+0.22 

1.6+0.210-'° 1 

2.6+0.5x10-'° 1.96+0.64 

i.5+0.5x10-'° 1 

1.9x10-° 4i) 

4.7x10-'° 1 


l-B(l) 


Transition 

(1-0) 

X^^ Moi 

/HCO+ 


HCN 

HCO+ 

1.3±0.2 10-‘“ 1.01+0.14 

1.3±0A 10-'° 1 



Notes. For each molecule, we provide the expected abundances, assuming LTE and optically thin conditions, to compare them 
with Galactic values. We provide abundances and abundance ratios for two positions M31-I(3) and M31-I(2) (columns 2-3) and 
the stacking detections (columns 4-5), and then compare these values with those expected for the Galactic interstellar medium, at 
4kpc (MR) and at the Galactic centre (GC), derived from Bergin et al. (1995) and Wilson & Rood (1994) (columns 6-9). When 
'^CO(l-O) is not measured directly, the numbers are in italics. 


which we adopt an ad-hoc estimate of based on the 
averaging of the measured and loco measurements 
performed with the spectra used for the stacking: a scaling 
factor is derived from these averaged measurements and 
applied to the stacked '^CO(2-l) measurement to estimate 
'^C(9(2-l). It is then scaled to '^CO(l-O) with the 0.8 fac¬ 
tor, as discussed in Sect. 3. For the sake of discussion, we 
also provide (in italics) M31I-B(1) and the stacking mea¬ 
surements for which we have HCN(l-O) and HCO4-(l-0) 
measurements to no direct '^CO(l-O) measurements. As 
shown in Table 8, the molecular hydrogen column den¬ 
sity exhibits variations up to a factor 3 owing to different 
resolutions (11" and 21"). 

We use the column densities presented in Table 4 and 
described in Sect. 3. We also provide for comparison the 
abundance values for the Galaxy computed in the molec¬ 
ular ring (MR) at 4kpc from the Galactic centre and in the 
Galactic centre (GC) by Bergin et al. (1995) and Wilson 
& Rood (1994). 

The abundances are affected by the uncertainties in 
the molecular hydrogen column densities computed em¬ 
pirically with the Xco factor and by the uncertainties in 
the excitation temperatures. For the CO transitions, we 
can expect that variations of the excitation temperatures 
are a secondary effect, and we can thus rely on abundance 
ratios for optically thin transitions. 

The detected ^^C0/C'*0 abundance ratios are very 
different from the Galactic values. There is a discrepancy 
of a factor 4 with the molecular ring (MR), and a fac¬ 
tor 9 with respect to the Galactic centre (GC) where gas 
has been processed regularly. Obviously, the gas is not 
in the same state as these obvious regions in the Galaxy. 


Considering the C^^O abundances in M31-I(3), it is rel¬ 
atively close to the molecular region value, while the 
discrepancy affects the '^CO. For the stacking, the C**0 
abundance is a bit smaller and intermediate between GC 
and MR values, even though it is compatible with the 
M311(3) abundance within 1 cr. While the C'^O abun¬ 
dance lies in the same range as the Galactic value, the 
'^CO abundance is deficient. 

The HCN and HCOh- abundances differ from the 
Galactic values, and the abundance ratio. These dis¬ 
crepancies could be due to different excitation temper¬ 
atures. However, our different measurements are glob¬ 
ally compatible within 3cr. Liu et al. (2013) measured 
an HCN/HCOh-( 1-0) abundance ratio of 1.2 + 0.4 in 
infrared dark clouds (IRDC) and discuss that this ra¬ 
tio may be unaffected by the environment. In parallel, 
they measured an average HNC/HCN abundance ratio 
of 1.47+0.50 in IRDC, while the 3cr upper limit we de¬ 
rived on the HNC/HCN abundance ratio for M31I-B(1) 
and M311(2) is 1.01 and 1.14. As these authors discuss, 
we can argue that our detections are different from the 
IRDC. Our upper limits allow us to put some constraints 
on the kinetic temperature (e.g. Schilke et al. 1992; Hirota 
et al. 1998). This corresponds to kinetic temperatures in 
the range 38-100 K and 14-50 K, and is compatible with 
the results of Viaene et al. (2014) based on the analysis of 
infrared data. Similarly, we get HCO+/HOC+ abundance 
ratios larger than 1.7 and 4.9, but this is not restrictive 
enough as previous detections are in the range 50-12500 
(Savage & Ziurys 2004). 

In the next section, we explore with non-LTE simu¬ 
lations how different excitation temperatures could affect 
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Fig. 11. Main physical parameters corresponding to the HCN gas observed in M31-I(2) and M31-I-B(l). The optical depths are 
shown as function of the HCN column density (per kms^*) and the molecular hydrogen density iiu^. The corresponding contours 
for these two clumps are shown with a dashed line for M31-I-B(l) and with a full line for M31-I(2). The vertical (green) line 
indicates the critical density computed for a collisional temperature of 20 K (see Table 3). The lower horizontal full (dashed) 
line shows the column density measured for M31-1(2) (M31-I-B(l)) at 9K. The upper horizontal lines correspond to the column 
densities computed with standard abundances (Bergin et al. 1995) and the corresponding molecular hydrogen column densities. For 
M31-I-B(l), we derive an upper limit as the molecular hydrogen column density is only available for a smaller beam (see Table 4) 
and we expect some dilution of the signal. The contour levels correspond to excitation temperatures of 9 K (full line), 15 K (dashed 
line), 20 K (dotted line), and 25 K (dash-dotted line). 


the different parameters and help to get a coherent set of 
parameters. 

4.5. RADEX simulations 

In order to check the consistency of our parameters, we 
ran some simulations with the radiative transfer code 
RADEX (van der Tak et al. 2007) for several molecules. 
The idea is to relax the LTE conditions to scan a range of 


excitation temperatures and of molecular hydrogen densi¬ 
ties. Knowing the column densities of each emission line, 
we can estimate the excitation temperature and the den¬ 
sity. As indicated in the figures discussed below, the stan¬ 
dard abundances correspond to optically thick conditions 
and can be excluded. 

In Eig. 11, we show the results of the simulations run 
for HCN and superimpose the various parameters dis¬ 
cussed in the previous section, namely, the optical depth 
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Fig. 12. Main physical parameters corresponding to the '^CO(2-l) gas (left panel) and the C'*0(2-l) gas (right panel) detected 
in M31-I(3). The optical depths are shown as function of the '^CO(2-l) column density (per kms^*) and the molecular hydrogen 
density «a 2 - The vertical green-black lines indicate the critical densities computed for a collisional temperature of 20 K (see Table 
3). The horizontal thick, full yellow-black lines indicate the '^CO(2-l) and the C**0(2-l) averaged column densities (per unit 
velocity), while the horizontal thin black lines indicate the corresponding 1 cr standard deviation; the hatched area corresponds to 
the M31-1(3) detection. The horizontal thick, dashed yellow-black lines (thin black lines) correspond to the column densities (1 cr 
standard deviation) computed with standard abundances (Bergin et al. 1995) and the corresponding molecular hydrogen column 
densities. The full line (dashed and dotted) contour levels correspond to an excitation temperature of 17 K (14 K and 19 K). The 
optical depths, corresponding to 17 K, are shown as a thick line contour level. 


and the molecule column densities per unit velocity de¬ 
rived from our observations and from standard abun¬ 
dances. These results support the view that this gas is sub- 
thermally excited at 9 K with a density smaller than the 
critical density, which is provided in Table 3. This critical 
density plotted here (and in Figs. A.5 and 12) does not 
match Tex - 2QK, confirming the arguments of Shirley 
(2015) that this quantity is an oversimplified and incorrect 
interpretation of the excitation density of a transition. The 
results of the constraints of RADEX simulations with our 
observations are also consistent with the HCN(l-O) tran¬ 
sition being optically thin. We detected two dense clumps 
with a molecular hydrogen density = 4.8 x 10^ cm“^ 
on the north-west side of the inner ring. Our abundances 
are more than a factor of 10 lower than the column den¬ 
sities derived from the standard abundances estimated by 
Bergin et al. (1995) in the molecular ring (MR), and the 
discrepancy would be even larger if one considers the 
Galactic centre conditions. 

We performed similar computations for HCOh-(I-O), 
as shown in Figure A.5. We also find subthermal con¬ 
ditions with a molecular hydrogen density smaller than 
the critical density and an excitation temperature of 9 K. 
The abundances are also smaller (factor 5-6) than the 
column densities derived from the standard MR abun¬ 
dances (Bergin et al. 1995). The gas associated with the 
HCO-i-(l-O) detection has a density of 1.0 x 10^ cm“^. 

We applied the same procedure to M31-1(3) where 
*^CO and C'^O have been detected, as shown in Fig. 12. 


On the one hand, the best excitation temperature is consis¬ 
tent with 17 K but we cannot exclude 20 K given the error 
bars. On the other hand, the simulations support a molecu¬ 
lar hydrogen density larger than the critical density. These 
lines are close to FTE conditions and trace clumps with a 
density larger than 2.24x10^ cm~^. The '^CO measured 
abundance is smaller than the standard abundances by a 
factor 5, while the C'^O measured abundance is com¬ 
patible within the error bars of the standard abundances 
estimated in the molecular ring. 

The abundances we have derived are systematically 
smaller than the standard abundances (Bergin et al. 1995) 
for the molecular ring, but for C'^O. This is consistent 
with our previous discussion that '^CO is optically thick 
and '^CO depleted, with FTE and optically thin condi¬ 
tions for and C'*0. HCN and HCO-i- abundances 
are consistent with subthermal excitation conditions with 
excitation temperatures of 9 K and are smaller than the 
standard abundances. As shown in Figs A.3 and 10, the 
HCO-i/HCN line ratio measured in the clouds studied 
here is similar to that measured in M31 ’s disc by Brouillet 
et al. (2005). 


5. Discussion 

In this section, we discuss our dense gas detection, how 
it can improve our knowledge of this inner region, and 
whether the properties of this gas are consistent with our 
understanding of previous observations. In Sect. 5.1, we 
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Fig. 13. Azimuthal (1” and 3''‘‘ columns) and spatial (2"'^ and 4''') distribution of X-ray sources (top) and planetary nebulae and UV 
sources (bottom). The top left panels show the deep Chandra-detected X-ray sources from Stiele et al. (2011). The top right panels 
show the X-ray sources from Hofmann et al. (2013). The bottom left panels indicate the planetary nebulae from Ciardullo et al. 
(1989) (see also Melchior & Combes (2013)). The bottom right panels show the HST UV sources from Rosenfield et al. (2012). 
The dashed line exhibits the best-fit profile from Kaaret (2002) corresponding to a disc with a position angle PA= 52 deg. In the 
right panels, the two minima are roughly at the same level, while in the left panels the first is deeper as indicated qualitatively with 
the dotted line. The ticks in the U' and 3''“^ figure columns correspond to the position angles of our detections. 


discuss the geometrical configuration of the gas: the in¬ 
formation derived from dust extinction affecting point 
sources, and the properties of various gas components 
observed at different wavelengths. In Sect. 5.2, we dis¬ 
cuss our knowledge of the heating sources and the proper¬ 
ties of the interstellar medium in this region. In Sect. 5.3, 
we discuss the kinematical information obtained with our 
millimeter observations. 

5.1. Geometrical configuration of the gas 

M31 has been targeted by different surveys. In Sect. 5.1.1, 
we study the distribution of point sources from different 
catalogues in search of possible signs of extinction that 
trace the gas. In Sect. 5.1.2, we summarise the arguments 
we can derive from the diffuse components. 

5.1.1. Point sources and extinction 

In Fig. 13, we considered four catalogues: X-ray sources 
from Stiele et al. (2011), deeper X-ray sources from 
Hofmann et al. (2013), planetary nebulae from Ciardullo 
et al. (1989), and UV sources from Rosenfield et al. 
(2012). The UV sources studied by Rosenfield et al. 
(2012) trace the old stellar population, and they have 
the same radial distribution as low-mass X-ray binaries, 
which correspond to the majority of the X-ray sources in 
Hofmann et al. (2013). The stellar kinematics has been 


studied by Saglia et al. (2010) with long slits: bulge stars 
exhibit a rotation with a maximum velocity of 100 km s“' 
along the major axis (considered at PA= 48 deg) and a 
zero (systemic) velocity along the minor axis). These ve¬ 
locities clearly do not correspond to the gas component 
we detected. 

For each catalogue, we plotted the azimuthal distri¬ 
bution of point sources within a radius of 300, 200, and 
140", depending on the area covered. These populations 
all belong to the same ellipsoid of PA^ 52 deg as first 
discussed by Kaaret (2002). The maxima of the sinusoids 
correspond to the major axis, while the minima (where 
our dense gas has been detected) correspond to the minor 
axis. 

The number of sources is significantly larger for 
Hofmann et al. (2013) than for Stiele et al. (2011); the 
former work is more sensitive as the data are stacked 
over a longer period. The work of Rosenfield et al. (2012) 
based on HST data contains an unprecedented number of 
sources thanks to deep integration and a gain in spatial 
resolution. We can observe that in these deep observations 
the two minima of the sinusoid are about at the same level, 
while the previous observations exhibit an asymmetry. 
This asymmetry is not really significant when compared 
to the Poisson error bars shown in the Fig. 13. While in ra¬ 
dio we do not expect much extinction, the UV sources in 
Rosenfield et al. (2012) should be very sensitive to extinc¬ 
tion. However, we do not detect any strong asymmetry: if 
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we compute the ratio of the numbers of sources from the 
near and far sides, we find 0.9676 (0.9662) if we con¬ 
sider 52 deg (77 deg). These ratios are significant as they 
have a Poisson noise of 0.012. If this asymmetry is due 
to extinction, it affects 4% of the studied area. This order 
of magnitude is consistent with the gas beam filling fac¬ 
tors computed in Sect. 4.3. However, we do not have the 
sensitivity to derive whether this patchy dust component 
belongs to the large scale disc (PA=77 deg) or the more 
face-on component detected in the stellar bulge popula¬ 
tion (PA= 52 deg). We can exclude the possibility that 
young stellar clusters from the disc component could cre¬ 
ate this inhomogeneity, as no HII region is detected in this 
area (Azimlu et al. 2011) beside the cluster close to the 
centre (Lauer et al. 2012). The (young) stellar UV cluster 
(Kang et al. 2012) closest to our detection has an age of 
1 Gyr and is at the offset (-111",152") from the centre 
and at a closest separation of 71 arcsec (271 pc in projec¬ 
tion) to our north-western dense gas complex. Dong et 
al. (2015) detected an intermediate age (300 Myr to 1 Gyr 
old) stellar population at solar metallicity associated with 
the bulge, which could contribute up to 50% of the V in¬ 
tensity. However, they only studied the radial dependency, 
while we expect a smooth spatial distribution as the stars 
have performed numerous galactic rotations since their 
birth. 

Assuming a uniform distribution in the stellar ellip¬ 
soid, one can estimate the ellipticity as done by Kaaret 
(2002). For the Rosenfield et al. (2012) (Hofmann et al. 
2013) sample studied here within 140" (300"), we esti¬ 
mate e - 0.24 (e = 0.44). This is consistent with the 
measurements of the radial variation of the ellipticity and 
position angle performed on J-2MASS photometry by 
Melchior & Combes (see Figure A.l 2013). This supports 
our previous finding that there is very little extinction. 

While these populations trace the bulge, there is an 
unexpected X-ray source alignment in the Hofmann et 
al. (figure 14, 2013) catalogue very close to the cen¬ 
tre. Similar to that observed in the outskirts of M31 (e.g. 
McConnachie et al. 2009), these X-ray sources could be 
part of a tidal stream whose stellar population could be 
hidden by the crowding at other wavelengths. The facts 
that it corresponds to an old stellar population and that 
it is aligned are strong arguments to argue that they are 
far from the centre and seen in projection. However, it is 
not aligned with our positions nor any diffuse components 
discussed below. 

5.1.2. Gas components 

While the stellar population mainly belongs to the bulge, 
the diffuse components, apart from the UV to near-IR stel¬ 
lar continuum, are gaseous. 

Radio component The analysis of the radio emission 
(GieBiibel & Beck 2014) reveals complementary infor¬ 
mation. There is a significant continuum flux at 3.6 cm 
and 6.2 cm, and a weak polarised component (~10%). 


There is a relatively flat spectral index while 

the contribution from the black hole is very weak and 
localised (Li et al. 2011). This flat spectral index can¬ 
not be accounted for by an inverted synchrotron emis¬ 
sion due to a present AGN, but we cannot exclude relic 
emission from past AGN activity (e.g. Ishibashi & Fabian 
2015). The flat spectral index could be due to the com¬ 
bination of synchrotron and free-free emissions. In the 
Block et al. (2006) scenario, the inner ring forms and ex¬ 
pands to its actual position in ~ 150 Myr. The correspond¬ 
ing past starburst activities (150-200 Myr old) could ac¬ 
count for the continuum emission stronger close to the 
centre. The polarisation map of GieBiibel & Beck (2014) 
is illustrated as contours in Fig. 14. The polarisation is 
blurred by the Faraday rotation of the electrons. On the 
north side, no polarisation is detected, while it is detected 
on the south side. This polarisation pattern corresponds to 
a good approximation to the simple modelling proposed 
in Melchior & Combes (see Fig. 16 and Tab. 7, 2011) 
of the inner ring with a PA of -35 deg and an inclina¬ 
tion of 40 deg, together with the two discs components 
strongly depolarising the north-western side. In addition 
in the south side, we clearly see some inhomogeneities in 
the polarised emission detected: there is a gap that cor¬ 
responds to the Ha(-i-[A//]) emission. This confirms the 
parameters of the inner ring’s location on the near side, 
i.e. in the south-east direction. This supports some pre¬ 
vious tentative interpretations (Ciardullo et al. 1988) that 
we are observing structures more face-on than the 77 deg 
main disc with the near side on the north side and a more 
edge-on structure on the south side. 

Supporting the post-starburst scenario, the central 
field exhibits a remmant tidal swirl observed in radio at 
3.6 cm and 6.2 cm (GieBiibel & Beck 2014) and in ionised 
gas (Ciardullo et al. 1988). These radio emissions could 
trace relic activities of stellar populations younger than 
200 Myr (Condon 1992) i.e. following the impact and ex¬ 
pansion of the inner ring. This remnant tidal swirl could 
be dominated by shocks, e.g. from supernovae as [Nil] is 
larger than H^ (Rubin et al. 1972), while the inner ring 
is possibly stronger in Hq, tracing a 200Myr old starburst 
(Block et al. 2006). 

Ionised gas While we know that the observed extinction 
is weak and the overall amount of gas is small, Melchior 
et al. (2000) has estimated a lower limit on A^ assum¬ 
ing all the gas is in front of the bulge. Even in loca¬ 
tions where the extinction detected is small (i.e. proba¬ 
bly underestimated by the previous assumption), it corre¬ 
sponds well with features detected in the infrared (e.g. at 
100jum as indicated in Figure 1). As discussed by Ford 
et al. (2013), the central gas is hotter than the inner ring. 
The lOOjum map exhibits hot dust and cold dust, while 
the 500 jjm map mainly shows cold dust in the inner ring 
(close to our observed positions). This central hot compo¬ 
nent corresponds to the ionised gas (Ciardullo et al. 1988). 
The ionsed gas is very inhomegeneous and its distribu¬ 
tion does not exhibit any clear sign of extinction. This is 
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compatible with the fact that the gas is very clumpy with 
a filling factor smaller than 1 %and that the near/far-side 
asymmetry, possibly due to extinction, is at the level of 
4%. 

Molecular gas Our dense gas detections discussed 
here exhibit a surprising velocity pattern. This is the 
wavelength range where the best velocity resolution is 
achieved (few km/s), and we observe in two positions 
along the minor axis where the systemic velocity is ex¬ 
pected atvelocities between -530 and -140 kms^'. This 
kinematical information demonstrates the presence of 
several molecular gas components along these lines of 
sight. In contrast to HI measurements (Braun et al. 2009; 
Chemin et al. 2009), we do not have a priori the sensi¬ 
tivity to detect molecular gas in the warp or the external 
parts of the galaxy. 

5.2. Heating sources 

We investigate here the possible heating sources, which 
could explain why the signal is stronger on one side 
and what heats the gas in general. Following Bogdan & 
Gilfanov (2008), we know that the diffuse X-ray gas de¬ 
tected in this area within 150-200" from the centre is 
weak. The detected 0.5-1.2keV flux lies in the range 
1.6 - 9xl0“'*erg s“'cm“^arcsec“^, considering a molecu¬ 
lar cloud of 1 pc with a distance of M31 of 780 kpc, we 
can estimate a rough X-ray incident radiation intensity 
of 6-34x10“^ ergs^'cm^^. This is at least five orders of 
magnitude weaker than the radiation field considered by 
Meijerink et al. (2007) for XDR modelling. There is only 
one ROSAT (Second ROSAT PSPC catalogue) source 
within 12"(22") from M31I (M31I-B) with a 14arcsec 
positional error, with a flux of 10“^^ erg s“^cm“^arcsec“^. 

There is an ultra-luminous X-ray source (Kaur et al. 
2012) located 22" from M31-2c-a, which has exhibited 
an outburst that lasted three months. ULX are corre¬ 
lated with star formation (Mapelli et al. 2009). As in the 
Cartwheel, this ULX is detected next to the inner ring. 
The presence of ULX in the ring supports the starburst 
scenario, as several ULX are detected in ring galaxies sus¬ 
pected of frontal collision (Prestwich et al. 2012). ULX 
are expected in low-metallicity environments (Prestwich 
et al. 2013), which could explain why only one has been 
detected. Barnard et al. (2012) measured an outburst flux 
of 50xl0^^erg s ', which could contribute at the level of 
5 X 10“^ erg s“' cm“^. 

In this inner region, the contribution of the cen¬ 
tral black hole could be significant. However, Li et al. 
(2011) measured a burst of activity in January 2006 with 
Chandra X-Ray Observatory at a level of 4.3x 10^^ erg s“' 
and a quiescent level of 4.8x10^^ ergs^'. At the dis¬ 
tance of M3II and M31I-B (supposed in the sky plane), 
this would contribute to the X-ray incident radiation in¬ 
tensity at a level of 6.9 and 5.5xl0“*erg s“'cm-2 (6.2 
and 5.0x10“^erg s“' cm“^) in the quiescent state (during 
the outburst). This is four to five orders of magnitude 


weaker than the possible contribution of the ULX dis¬ 
cussed above. 

The X-ray sources shown in Fig. 14 are strongly con¬ 
centrated towards the centre. The catalogue Hofmann et 
al. (2013) based on Chandra HRC is more sensitive than 
other catalogues (Stiele et al. 2011; Barnard et al. 2014). 
The surface density varies as (R/1 kpc)“^ ^. The contribu¬ 
tion of the central part is thus more important than the 
black hole activity. 

We can thus exclude that the excitation mechanism of 
the gas is due to X-ray radiation from the gas accreted 
by the black hole. The [OIII] line emissions detected by 
Saglia et al. (2010) and M. Sarzi (private communication) 
can be accounted for by planetary nebulae and no LINER 
activity is necessary. 

The heating is mainly due to the old stellar popula¬ 
tion as discussed by Draine et al. (2014) and Viaene et al. 
(2014). Transient X-ray sources, such as ULX, could con¬ 
tribute inhomogeneously to the heating. The heating rate 
U, based on Mathis et al. (1983) and computed by Draine 
et al. (2014), are provided in Table 1. The area around 
M31I and M31I-B have U — 1.5, while the other three 
positions have U in the range 1.8-2.7. This corresponds to 
an average dust temperature of 20 K, according to Draine 
et al. (2014). 

5.3. Different veiocity components 

The kinematics of this area is complex with up to three 
different velocity components. The probability that the 
warp component, coming from radii larger than 15- 
20 kpc, is seen in CO emission towards the centre is neg¬ 
ligible. This is the case because of the weak metallicity 
of this gas, which makes it undetectable in CO. As dis¬ 
cussed in Melchior & Combes (2011), if we consider that 
the probability of chance alignment of gas is low and that 
the gas is inside the main galaxy, this means that we are 
in the presence of several structures. Indeed, an explana¬ 
tion in terms of non-circular motions in a possible barred 
potential has already been ruled out given the velocity 
disposition and the fact that the gas distribution is not 
aligned with any elongated structure, but in an off-centred 
ring (Melchior & Combes 2011). The component at the 
systemic velocity is compatible with gas along the mi¬ 
nor axis of the main disc, while the two counter-rotating 
components could be orbiting in two planes with differ¬ 
ent inclinations and orientations. They could correspond 
to an inclined disc and to the inner 1-kpc ring. However, 
in presence of observations in two positions, it is diffi¬ 
cult to fully constrain such a configuration. Apart from 
the inner ring detected in infrared (see Fig. 1 and Block et 
al. (2006)) and the Hq,(-i-[A//]) structure detected inside 
(Ciardullo et al. 1988), we do not have kinematic obser¬ 
vations covering these positions. Boulesteix et al. (1987) 
got the [Nil] velocity field within a 3'x3' field and Nieten 
et al. (2006) did not integrate deep enough in the central 
part to detect CO. Opitsch et al. (2015) has detected two 
rotating components in [OIII], but their two components 
rotate in the same direction, which adds complexity to the 
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Fig. 14. Known heating sources. On the left side, we superimpose on the Htr+iNII] map from Ciardullo et al. (1988): the X-ray 
sources detected by Hofmann et al. (2013) with blue dots; the 1 Gyr young stellar cluster detected by Kang et al. (2012) with the 
green triangle; and the X-ray sources detected by Stiele et al. (2011) and Barnard et al. (2014) with white and purple dots. The 
blue dots with a yellow circle correspond to the first (and more central) ultra-luminous X-ray transient (ULX) detected in M31 (e.g. 
Kaur et al. 2012). The black contours correspond to the apparent As extinction map discussed in Melchior et al. (2000). The light 
blue contours correspond to the 6 cm polarised intensity detected by GieBiibel & Beck (2014). The large red circles correspond to 
our detections. Light black circles provide typical projected radial distance from the centre (x pc and y pc). Thick ellipses show the 
approximate position of the inner ring. On the right side, we indicate, in the top panel, the number of X-ray sources (Hofmann et 
al. 2013) detected at each radius; in the bottom panel, the X-ray luminosity at each radius is shown as a function of the radius. The 
mean and median values of the X-ray luminosity (averaged in lOOpc radial bin) is also provided. 


problem. Interestingly, the polarisation map of GieBiibel 
& Beck (2014) corresponds to the parameters proposed 
by Melchior & Combes (2011) for the three-component 
modelling. 

Finally, the amplitudes between our different veloc¬ 
ity components are different and larger than the simulated 
modulation proposed by Dierickx et al. (2014). 

The alternative interpretation is that we are observing 
gas superimposed on the central part of M31, i.e. gas in 
the outskirts of the disc (Koch et al. 2015; Lehner et al. 
2015; Wolfe et al. 2013; Lockman et al. 2012) is ruled 
out. No high velocity clouds have been detected in CO 
(Combes & Charmandaris 2000; Miville-Deschenes et al. 
2005), while Galactic gas should be around 0kms '. 

5.4. C depletion and evidence for a post-starburst 
scenario 

In Sect. 4.1, we discussed that it is remarkable that the 
C'^O detection is at a level comparable to the '^CO de¬ 
tection, in contrast to what is observed in the Galaxy (e.g. 
Nishimura et al. 2015; Shimajiri et al. 2014), where the 
C'^O is usually 10 to 150 times weaker than CO. In 
Table 7, we have shown that '^CO and C'^O are optically 
thin, so we can claim that the *^CO is under-abundant 
with respect to C'^O. In Sect. 4.5, we have shown that 
these two CO components are close to thermal equilib¬ 
rium. Hence, we can rely on the abundance ratios com¬ 
puted in Table 8; the deficit is five to ten times weaker 
than with the Galactic abundances. This deficit could be 


accounted for by a post-starburst scenario. This kind of 
*^CO deficit was first observed by Casoli et al. (1991) in 
post-merger. Davis (2014) also find an anti-correlation of 
*^CO and star formation and gas surface densities in local 
galaxies. 

This could be the signature of a 200 Myr old star- 
burst triggered by a collision with a galaxy like M32, as 
proposed by Block et al. (2006). This scenario accounts 
for the two off-centred ring structures of M31, with most 
star formation in the lOkpc ring (Ford et al. 2013) and a 
smaller ring, where we have detected CO and dense gas. 

6. Conclusion 

We have detected *^CO(2-l) and *^CO(2-l) lines with 
multiple velocity components along the minor axis of 
M31 on both sides of the inner 1-kpc ring. We detected 
dense gas traced by C'*0(2-l), HCN(l-O), and HCO-h( 1- 
0) on the north-west side at the systemic velocity for 
C'®0(2-1) and in the blue-shifted component for HCN(1- 
0) and HCOh-(I-O). We also detected a redshifted compo¬ 
nent for the three molecules by stacking at the CO veloc¬ 
ity on both sides. 

Draine et al. (2014) and Viaene et al. (2014) found 
lower measurements for the stellar heating in this area 
on the north-west side, which could explain why dense 
gas has only been detected directly on this side. A careful 
analysis has shown that '^CO is optically thick while the 
other CO lines are optically thin and close to LTE at an ex¬ 
citation temperature of 17.5-20 K. The gas is very clumpy 
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and the averaged beam filling factor is 0.8%. The average 
column density of molecular hydrogen is 16 lO^^cm^^. 
The derived abundances are all smaller than the stan¬ 
dard abundances, but for C'*0. The comparison of the 
derived abundances with RADEX simulations shows that 
the HCN and HCO-h gas correspond to dense clumps with 
molecular hydrogen density of the range 1 - 6 x 10^ cm“^ 
and an excitation temperature of 9 K. The HCO-I-/HCN 
line ratio is comparable with that measured for clouds 
detected in M31’s disc (Brouillet et al. 2005), while the 
derived abundances are smaller than the Galactic abun¬ 
dances (Bergin et al. 1995), with a larger discrepancy for 
HCN(l-O). Only one transient ULX object has been iden¬ 
tified as a possible secondary source of heating, while the 
bulge stellar population is sufficient to heat the gas and 
reproduce the infrared maps. 

Our isotopic CO lines and HCN/HCOh- lines detec¬ 
tions are all weak and close to the detection limits. We 
detected C'^O in a single point at the systemic velocity. 
Stacking of all '^CO detections has revealed the under¬ 
lying presence of C'^O gas in the redshifted component. 
Both the single point and the stacking detections have an 
overall intensity comparable to '^CO. As we expect C'*0 
to be six times less abundant than ^^CO and the measured 
abundance of C^^O is compatible with the standard abun¬ 
dance, ^^CO should be depleted. A coherent explanation 
is that has been enriched during the recent staburst 
and ^^CO had no time to catch up, so it appears depleted. 
This is compatible with the Block et al. (2006) scenario 
suggesting a frontal collision with a nearby companion 
galaxy such as M 32. This collision explains the rings and 
has triggered a starburst at the epoch of the collision. This 
is also in agreement with the central starburst, which oc- 
cured next to the black hole 200 Myr ago (Lauer et al. 
2012). 

We confirm the presence of three kinematical com¬ 
ponents detected in molecular gas, as first discussed in 
Melchior & Combes (2011). Separated by 400km s“^ and 
on both sides of the systemic velocity, this configura¬ 
tion on the minor axis is not expected for a disc in ro¬ 
tation and cannot be accounted for by a bar. We have also 
checked that the parameters we first proposed for a simple 
modelling of the three kinematical components are com¬ 
patible with other observations. In addition to the 77deg 
main disc with a position angle (PA= 35 deg) detected 
in HI, there is an less inclined nuclear disc (i= 43 deg) 
with PA= 53 deg (rather than 70 deg, as initially thought). 
The position angle of this nuclear disc is well defined 
with different stellar sources catalogues (Kaaret 2002; 
Rosenfield et al. 2012). The 1-kpc ring with i= 40 deg and 
PA= -35 deg accounts for the polarised map of GieBiibel 
& Beck (2014) very well. 
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Fig.A.2. HCN(l-O) and HCO-)-(l-O) stacking of the south¬ 
east positions. The stacking is based on the *^CO velocities. 
The top (bottom) panel shows the stacking corresponding to 
> -300km s“' (< -300kms“*) velocities. The velocities are 
shown in the rest-frame frequency of the HCO-l-(l-O) line. 
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Fig. A.l. Averaging of the observations performed in the EMIR lower outer (LO) bands. The different lines present in the bandwidth 
are plotted in the rest-frame frequency. No signal is detected. This has been used to derived detection limit in Table 4. 
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Fig. A.3. Comparison with other catalogues (2). HCO+/CO(1-0) line ratio versus HCN/CO line ratio. The two detections presented 
here (M3II and M31I-B) are superimposed on previous measurements, namely: detections across M3rs main disc by Brouillet et 
al. (2005) and a compilation of various types of galaxies from Krips et al. (2010). The *^CO(2-l) measurements presented here have 
been corrected by a CO 2-1/1-0 line ratio of 0.8, as measured in Melchior & Combes (2011) in a complex M31G within 10 and 30 " 
from the observed positions. The line indicates the best fit provided by Brouillet et al. (2005) for M31’s main disc. Our points are 
systematically above the best fit they provide, suggesting an excess of HCO+. The overall distribution of Krips et al. (2010) does 
not exhibit any clear pattern. 
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Fig.A.4. C'^O/'^CO (2-1) line ratios versus the '^CO/*^CO(2-l) line ratios presented in this paper. Arrows indicated 3(T upper 
limits. The full black (empty yellow) point with error bars corresponds to our direct detection (our red stack detection). The line 
indicates a slope 1. 


25 














Anne-Laure Melchior, Franjoise Combes: Dense gas tracing the collisional past of Andromeda 


Hcrit Log,0 T 



3.75 4 4.25 4.5 4.75 5 5.25 5.5 5.75 6 

Logjo Hjj 2 


Fig.A.5. Main physical parameters corresponding to the HCO+(1-0) gas observed in M31-1(2) and M3 l-I-B(l). The optical depths 
are shown as function of the HCO-l-(l-O) column density (per km s“‘) and the molecular hydrogen density 11^2 ■ The vertical (green) 
line indicates the critical density computed for a collisional temperature of 20 K (see Table 3). The lower horizontal full (dashed) 
line shows column densities measured for M31-1(2) (M31-I-B(l)) at 8K. The upper horizontal line corresponds to the column 
densities computed with standard abundances (Bergin et al. 1995) and the corresponding molecular hydrogen column densities. For 
M31-I-B(l), we derive an upper limit as the molecular hydrogen column density is only available for a smaller beam (see Table 4) 
and we expect some dilution of the signal, as is the case for M31-1(2). The contour levels correspond to excitation temperatures of 
9 K (full line), 10 K (dashed line), 15 K (dotted line), and 20 K (dash-dotted line). 
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